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Putative global minima have been located for clusters interacting with an aluminium glue potential 
for N < 190. Virtually all the clusters have polytetrahedral structures, which for larger sizes involve 
an ordered array of disclinations that are similar to those in the Z, H and a Frank-Kasper phases. 
Comparisons of sequences of larger clusters suggest that the majority of the global minima will 
adopt the bulk face-centred-cubic structure beyond A'^ « 500. 
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I. INTRODUCTION 

The understanding of the structure of metal clusters 
has seen many developments in recent years/'^ however 
there is still much to be learnt. From the theoretical per- 
spective, for example, only relatively recently have global 
optimization techniques become sufficiently powerful to 
find the most stable structures of metal clusters with up 
to 100 atoms, ^ even when described by relatively simple 
many-body potentials. This has led to many interesting 
new structures being revealed''"^^ that go beyond those 
classified for pair potentials, but there are still probably 
many classes of structure that remain undiscovered. Fur- 
thermore, the many-body nature of the bonding makes it 
difficult to unravel the relationship between the observed 
structures and the form of the interactions, although 
there has been some interesting recent progress. ^'^^'^"^ 
This situation contrasts with that for pair potentials, 
where, for example, the effects of the width of the po- 
tential well^"*'^^ and oscillations in the potential^^"^^ on 
the cluster structure have been mapped out and ratio- 
nalized. 

One of the more interesting recent discoveries has 
been that some monatomic metal clusters'*'"'^^ can ex- 
hibit structures with polytetrahedral order, i.e. the 
whole of the cluster can be naturally divided up into 
tetrahedra with atoms at their vertices. In the exper- 
iments of Dassenoy et al. on cobalt clusters, the elec- 
tron microscopy images and diffraction patterns could 
not be explained by any of the usual structural forms 
for a cluster, but instead were consistent with a polyte- 
trahedral structure. Furthermore, the magic number at 
N — 61 in the mass spectrum of strontium clusters^'^^ 
was assigned to a polytetrahedral structure,* and the 
next magic number at iV = 82 also coincides with a poly- 
tetrahedral magic number. Similar magic number pat- 
terns have also been found for some rare-earth metals.^* 

These results seem somewhat surprising since ordered 
polytetrahedral structures are restricted to alloys in 
bulk, e.g the crystalline Frank-Kasper phases^^"^^ and 
quasicrystals,^^ structural models of which often invoke 
substantial polytetrahedral order. However, as we 
will see, there are good reasons why polytetrahedral 
structures are more likely to be seen for monatomic sys- 
tems in clusters. 

Figure 1(a) illustrates that all space cannot be packed 




FIG. 1: (a) Packings of regular tetrahedra to illustrate the 
strain inherent in polytetrahedral packings. Five regular 
tetrahedra around a common edge produce a gap of 7.36°, 
and twenty regular tetrahedra about a common vertex pro- 
duce gaps equivalent to a solid angle of 1.54 steradians. (b) 
The Frank-Kasper polyhedra for coordination numbers 12, 14, 
15 and 16, as labelled. The associated disclination networks 
are placed next to each cluster. 



with regular tetrahedra. (For this reason polytetrahedral 
structures are sometimes said to be frustrated. ^^) There- 
fore, the nearest-neighbour pair distances in a polyte- 
trahedral structure cannot all take the same value. For 
example to close the gaps associated with the icosahe- 
dral packing of regular tetrahedra, the tetrahedra must 
be distorted such that the edges of the tetrahedra on 
the surface are 5% longer than the radial edges. For a 
monatomic system this strain usually leads to an ener- 
getic penalty, although for a binary system, if a smaller 
atom of the right size is placed in the centre of the icosa- 
hedron, all pair distances could still take their ideal val- 
ues. 

For clusters of sufficiently small size the energetic 
penalties associated with the strain inherent to polyte- 
trahedral structures can often be more than compensated 
by a favourable surface energy. The surface of a polyte- 
trahedral cluster is triangulated leading to high coordina- 
tion numbers, whereas clusters based on a close-packing 
scheme may have to expose {100} faces, which involve 
atoms with lower coordination number, and hence typi- 
cally have a higher surface energy. 

The 13-atom icosahedral structure, which is very com- 
mon for systems interacting with both pair potentials 
and many-body metal potentials, illustrates this prin- 
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ciple. The polytctrahedral packing of the icosahedron 
can be continued by adding capping atoms above the 
faces and vertices. This "anti-Mackay" overlayer leads 
to polyicosahedral structures, where each atom in the in- 
terior of the cluster has local icosahedral coordination. 
So, the 19-atom double icosahedron is also common, but 
as the size increases these structures become less likely 
to be observed. For example, for Lennard-Jones clusters 
they are lowest in energy up to — 30. ^'^ The reason 
that these clusters become disfavoured is clear from Fig. 
1(a). As further regular tetrahedra are added to these 
packings the gaps involved become larger and larger, and 
hence the strain energy associated with polyicosahedral 
packings increases very rapidly with size. 

For the above polytctrahedral clusters, there are five 
tetrahedra around each edge in the interior of the clus- 
ter. To balance out the increasing gaps in this type of 
packing, the bulk Frank-Kasper phases also involve sites 
where six tetrahedra share a common edge — a negative 
diselination is said to run along this edge.^'' Although 
the local distortion required to remove the overlap when 
six regular tetrahedra are placed around a common edge 
is larger, the introduction of diselination lines leads to a 
decrease in the overall strain. 

The coordination polyhedra involved in the Frank- 
Kasper phases^"^'^^ are depicted in Figure 1(b). Local 
icosahedral coordination involves no disclinations, but as 
the coordination number increases negative disclinations 
must be introduced. Note a diselination line can only 
end by forming a loop or exiting at a surface. So for 
the Z14 coordination polyhedron the diselination passes 
through the central atom, and for the Z15 and Z16 co- 
ordination polyhedra the central atom acts as a node for 
three and four disclinations, respectively. Polytetrahe- 
dral packings can therefore be described by a network of 
diselination lines threading an icosahedrally-coordinated 
medium. In Frank-Kasper crystals this diselination net- 
work is ordered and periodic, and is mediated by the 
Z14, Z15 and Z16 Frank-Kasper coordination polyhedra. 

Polytctrahedral clusters that involve disclinations have 
previously been found for two very different pair poten- 
tials. Firstly, for clusters interacting with a sufficiently 
long-ranged Morse potential the strain involved in a poly- 
tctrahedral structure can be tolerated, and ordered poly- 
tctrahedral clusters are found up to A'' ps 70.^^ Secondly, 
a modified Dzugutov potential favours clusters that are 
based on the C14 and C15 Frank-Kasper phases up to 
at least TV = 250.1*^ This potential has a maximum m 
the potential that has been chosen to coincide with the 
next-nearest neighbour distance in close-packed struc- 
tures, and somewhat resembles the Fricdel oscillations 
that can occur in effective pair potentials derived for 
metallic systems. ^^'^^ 

These two examples arc for somewhat unusual poten- 
tials. However, polytctrahedral structures are perhaps 
more likely for metals, because the nature of the metal- 
lic potentials makes them potentially less sensitive to the 
strain in the individual pair distances, for reasons that 



will be discussed in more detail in Section II A. For ex- 
ample, similar structures to those of the Morse clusters 
have also been found for clusters described by a gener- 
alized set of metallic interactions."^* Furthermore, in this 
paper I report that the lowest-energy structures of clus- 
ters interacting with an aluminium glue potential^® are 
also polytctrahedral up to at least N = 190. 

The paper is organized as follows. In Section II the 
main focus is on the form of the aluminium potential 
and the implications of this form for cluster structure. 
In Section III I describe in detail the structures of the 
putative global minima, and in Section IV I take Algi as 
a case study for understanding the reasons for this sys- 
tem's preference for polytctrahedral structures. Then, in 
Section V I briefly examine the size evolution of the clus- 
ter structure at larger sizes. The emphasis in the paper is 
less on the clusters as a realistic model of aluminium clus- 
ters that will hold up to detailed scrutiny — this is very 
difficult for any empirical potential — and more on using 
these clusters as a model system for exploring the types 
of polytctrahedral structure that metal clusters might ex- 
hibit, and the reasons for their formation. However, in 
Section V I do compare the results with the available 
experimental information on the structure of aluminium 
clusters. 



II. METHODS 

A. Potential 

To model the aluminium clusters I use a glue 
potential^" that has been constructed using the force- 
matching method. The potential energy is given by 

E = -Bpair + Egi-ae (1) 

i<j i 

where ^(r) is a short-ranged pair potential, U{p) is a 
many-body glue function and pi is defined as 

= E^ 

j 

where p{r) is an "atomic density" function. These three 

functions have been fitted to match the forces produced 
by first-principles electronic structure calculations for a 
large set of configurations, including those that corre- 
spond to surfaces, clusters, liquids and crystals. That 
the potential has been designed to model a wide vari- 
ety of atomic environments, such as occur in clusters, 
is particularly important for the current application, as 
opposed to a potential that has just been fitted to bulk 
data. This flexibility has also led the potential to be 
used in a wide variety of applications, often with con- 
siderable success. Examples include clusters, wires, ''^ 
surface thermodynamics,''"^'"'^ grain boundaries,''^ bulk''^ 
and surface^'' diffusion and melt solidification.''* 
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FIG. 2: The three functions that make up the aluminium glue 
potential: (a) <f){r), p{r) and (b) U{p). 



The functions ^(r), p{r) and U{p) are displayed in Fig- 
ure 2. The forms of these functions are not assumed, but 
are determined by the fitting process. This contrasts with 
other empirical potentials, such as the Gupta or Sutton- 
Chcn potentials. With these latter types of potentials 
it would be impossible, for example, to reproduce the 
double well structure of ^!>(r), which also occurs for a so- 
phisticated 'bulk' potential.''^ 

The shape of these functions have some fairly straight- 
forward consequences for the cluster structure. U{p) has 
a single minimum of depth —2.6571 cV at p = 1.0. This 
term provides the dominant contribution to the energy, 
so it is important that each atom has a value of p as 
close as possible to one. Of course, the values associ- 
ated with the surface atoms will be significantly less than 
this. There is therefore a tendency for surface distances 
to contract so that the atoms achieve as high a value of 
p as is possible. This effect will be balanced by the ener- 
getic penalty associated with the resulting compression 
of the interior and the reduction in the pair energy of the 
surface atoms. Furthermore, the minimum of the glue 
function is rather broad, and so the energy loss associ- 
ated with p values that deviate from one is not so severe. 



Thus, surface contraction will not be as pronounced as 
for some other metallic systems. For example, the lead 
glue potential,''" whose cluster structures have recently 
been analysed,''"''^ has a stronger dependence of U on p. 

p{r) decreases roughly linearly with distance until r « 
4 A, after which it decays more slowly to zero at the cutoff 
distance. As the value of p{r) in the latter regime is 
always small, the major contribution to the p values of 
each atom comes from nearest neighbours, and hence the 
average value of pi will be maximized by clusters with a 
large number of nearest neighbours. 

(j){r) has an interesting double-well structure. The first 
minimum occurs at rmini = 3.095 A, and has a depth of 
— 0.0587 eV. This is separated by a maximum at 3.669 A 
from a second minimum at r^in2 = 4.303 A = 1.390 rmini 
of slightly greater depth, — 0.0633 eV. As the number 
of next-nearest neighbours will be significantly larger 
than nearest neighbours (except for the smallest clus- 
ters), the largest contribution to impair comes from the 
second neighbour shell. 

The ratio of the positions of the two minima would ini- 
tially suggest that they are almost ideally placed to coin- 
cide with the first and second shells for closed-packed 
structures, which have a ratio of v^. However, the 
nearest-neighbour distance is determined more by the 
glue part of the potential. For example, for an atom with 
twelve nearest neighbours, these neighbours would need 
to be at a distance for 2.611 A to achieve a p value of 1. If 
the effects of next neighbours are also included, this dis- 
tance does not not need to be as short, and so the nearest 
neighbour distance for the equilibrium facc-ccntred-cubic 
(fee) crystal, r^^^, is 2.851 A.^^ Compared to this dis- 
tance the position of the second minimum is not so advan- 
tageous for closed-packed structures; rmin2 = 1.509 rj^jf'. 

One of the important features of the glue energy for the 
current study is that it does not directly depend on the 
distribution of pair distances or the coordination number 
of an atom, but only on the pi values. This situation 
contrasts with that for a pair potential, and makes metal- 
lic systems potentially less sensitive to internal strains. 
This feature has already been cited as an important fac- 
tor in the stability of disordered clusters for a variety 
of metallic systems. The same principle applies for or- 
dered polytetrahedral clusters, which have both strain 
and high coordination numbers. For example, the ideal 
glue energy can be achieved for an atom with coordina- 
tion number higher than twelve simply by expanding the 
nearest-neighbour distances for this atom with respect 
to rjjjf' until p = 1 for that atom. Similarly, there is no 
reason why an atom with a wide distribution of nearest- 
neighbour distances cannot also achieve a p value of one. 



B. Global Optimization 

The global optimization of the aluminium clusters was 

performed using the basin-hopping'^'"'^^ (or Monte Carlo 
minimization^"^) approach. This method has proved par- 
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ticularly successful in locating putative global minima 
for a wide variety of cluster systems. As the results 
are extended up to 190 atoms, a considerable compu- 
tational effort was required. It proved particularly im- 
portant to supplement the standard unbiased runs from 
random starting points for each size, with runs started 
from low-energy minima of nearby sizes with the appro- 
priate number of atoms added or removed. These latter 
runs were applied iteratively until no further new global 
minima were located. 

It should be noted that, of course, there is no guar- 
antee that I have been able to locate the true global 
minima, and the probability that a global minimum has 
been missed will increase with cluster size, as the size of 
the search space, and hence the number of minima,®^"^^ 
increases exponentially with N. However, from examina- 
tion of the statistics of how often independent runs locate 
the same lowest-energy minimum, I am confident that 
virtually all of the putative global minima up to iV = 100 
cannot be bettered, and that beyond this size the puta- 
tive global minima, if not truly global, are very good esti- 
mates. In this latter size range the runs seeded from low- 
energy clusters of other sizes become more important. 
Therefore, global minima are most likely to be missed 
when they are structurally different from those of nearby 
sizes. For example, the lowest-energy structure found in 
a basin-hopping run for AI147 was at —436. 6892 eV, even 
though the 147-atom Mackay icosahedron has an energy 
that is 0.0133 eV lower. However, this is the only exam- 
ple I found where reoptimization of a structure that is 
known to often be particularly stable bettered the global 
optimization results. I also tried to construct large poly- 
tetrahedral structures based upon some of the structural 
motifs present for N < 100, but again this did not lead 



to any improvements. 

III. GLOBAL MINIMA FOR N < 190 

The energies and point groups for the putative global 
minima are given in Table I. Point files are available 

online at the Cambridge Cluster Database."''^ The ener- 
gies of the global minima are represented in Figure 3 in 
such a way that makes particularly stable clusters stand 
out. A selection of clusters are depicted in Figure 4 that 
are either particularly stable or have some interesting 
structural feature. Most of the intermediate sizes involve 
structures similar to those in this figure, but normally 
with a somewhat incomplete surface structure. 

The only other study that has examined the structures 
of aluminium clusters interacting with the present poten- 
tial studied clusters with 13, 43, 55 and 147 atoms. For 
A^=13, 55 and 147, Sun and Gong considered high sym- 
metry structures and found the Mackay icosahedra to be 
lowest in energy, in agreement with the results here. For 
A'^=43, they performed simulated annealing, generating a 
lowest-energy structure that is 0.0562 eV above the global 
minimum, and which is in fact the second lowest-energy 
isomer. They also claimed that this cluster had a 'glassy' 
structure, but in common with the global minimum it in 
fact has an ordered polytetrahedral structure. 

The structures with from 3 to 13 atoms are those 
typically seen for isotropic potentials. Most are on a 
polytetrahedral growth sequence leading to the 13-atom 
icosahedron, except for the six-atom octahedron and the 
8-atom dodecehedron. Beyond A'' = 13 growth occurs 
around the icosahedron with atoms being added above 
the faces and vertices of this structure in the so-called 
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TABLE I: Energies and point groups (PG) of the putative global minima. The unit of energy is eV. 



N 


PG 


Energy 


iV 


PG 


Energy 


N 


PG 


Energy 


N 


PG 


Energy 


N 


PG 


Energy 


3 




-4.099029 


41 


C2-0 


-112.178113 


79 


Cs 


-227.336269 


117 


Ci 


-343.817038 


155 


Cs 


-461.506096 


4 


Td 


-6.242292 


42 


Ci 


-115.154425 


80 


Cs 


-230.318346 


118 


Ci 


-346.764630 


156 


C-j,^ 


-464.682681 


5 




-8.607257 


43 


Cs 


-118.128134 


81 


Cs 


-233.343636 


119 


Ci 


-349.862568 


157 


Ci 


-467.796799 


6 


Oh 


-11.257920 


44 


C2v 


-121.260703 


82 


C2 


-236.308999 


120 


C2 


-353.034842 


158 


C2 


-470.980712 


7 




-13.860178 


45 


Ci 


-124.074147 


83 


Ci 


-239.278725 


121 


Ci 


-356.055763 


159 


Cs 


-474.154929 


8 


C2v 


-16.353276 


46 


Cs 


-127.238882 


84 


Cs 


-242.456890 


122 


C2 


-359.102294 


160 




-477.328706 


9 


Cl-a 


-19.063671 


47 


Cs 


-130.228220 


85 


C2v 


-245.646464 


123 


Ci 


-362.088084 


161 


Ci 


-480.408455 


10 


Czv 


-21.862051 


48 


Czv 


-133.418980 


86 


Cs 


-248.606413 


124 


Ci 


-365.266998 


162 


Cs 


-483.573938 


11 


C2v 


-24.616176 


49 


Cs 


-136.111570 


87 


Ci 


-251.592963 


125 


Ci 


-368.426994 


163 


Cs 


-486.653798 


12 




-27.780117 


50 


Cs 


-139.090832 


88 


Cs 


-254.736772 


126 


Ci 


-371.572700 


164 


C2 


-489.818114 


13 


h 


-31.278787 


51 


Ci 


-142.098016 


89 


Cs 


-257.685396 


127 


Ci 


-374.661505 


165 


Cs 


-492.896801 


14 




-33.585594 


52 


Ci 


-145.091286 


90 


Cs 


-260.675290 


128 


Ci 


-377.808061 


166 


Cs 


-496.054542 


15 


C2v 


-36.321872 


53 


Cs 


-148.261944 


91 


D3 


-263.877164 


129 


Cs 


-380.862378 


167 


Ci 


-499.123863 


16 


Cs 


-39.039888 


54 


Csv 


-151.376943 


92 


C2 


-266.936839 


130 


Ci 


-383.890586 


168 




-502.260998 


17 


C2 


-41.750455 


55 


h 


-154.612749 


93 


Ci 


-270.089293 


131 


Ci 


-387.030763 


169 


Ci 


-505.328765 


18 


Cs 


-44.777004 


56 


C-iv 


-157.245282 


94 


C2v 


-273.244125 


132 


C2 


-390.138328 


170 


Cs 


-508.459278 


19 


Dsn 


-48.182587 


57 


Si 


-160.381251 


95 


Cs 


-276.204093 


133 


Cs 


-393.177026 


171 


Cs 


-511.527235 


20 


C2v 


-50.823659 


58 


Cs 


-163.479187 


96 


C2 


-279.154395 


134 


C2 


-396.274720 


172 




-514.656897 


21 


Ci 


-53.470949 


59 


C2v 


-166.601625 


97 


Ci 


-282.115807 


135 


Cs 


-399.361340 


173 


C2 


-517.723146 


22 




-56.404346 


60 


Csv 


-169.702701 


98 


Ci 


-285.257294 


136 


Cs 


-402.433300 


174 


Ci 


-520.801927 


23 


Dsh 


-59.732308 


61 


Td 


-172.787060 


99 


Ci 


-288.206364 


137 


Cs 


-405.526322 


175 


Ci 


-523.966150 


24 


C2v 


-62.350544 


62 


Cs 


-175.512306 


100 


Ci 


-291.309715 


138 


C2v 


-408.694151 


176 


Ci 


-527.044129 


25 


Cs 


-65.195739 


63 


Ci 


-178.554122 


101 


Ci 


-294.334721 


139 


Cs 


-411.773939 


177 


C2 


-530.208354 


26 


Td 


-68.459809 


64 


C2v 


-181.598713 


102 


Ci 


-297.418136 


140 


Cs 


-414.942360 


178 


Ci 


-533.271452 


27 


C2v 


-71.162534 


65 


Ci 


-184.466611 


103 


Ci 


-300.517946 


141 


Cs 


-418.022390 


179 


Ci 


-536.411855 


28 


Cs 


-73.946989 


66 


Ci 


-187.498248 


104 


Ci 


-303.576439 


142 


C2v 


-421.159084 


180 


Ci 


-539.475366 


29 


D3h 


-77.115629 


67 


Ci 


-190.610516 


105 


Cs 


-306.773455 


143 


Cs 


-424.234289 


181 


Cs 


-542.615132 


30 


C2v 


-79.792682 


68 


Ci 


-193.508419 


106 


Ci 


-309.860784 


144 


C2v 


-427.293974 


182 


Ci 


-545.590132 


31 


Cs 


-82.634873 


69 


Ci 


-196.696161 


107 


Cs 


-312.991830 


145 


Ci 


-430.398146 


183 




-548.728221 


32 


C2v 


-85.678623 


70 


C2 


-199.701093 


108 


Cs 


-315.964908 


146 


Ci 


-433.528391 


184 


C2 


-551.764208 


33 


Csv 


-88.453125 


71 


Cs 


-202.700030 


109 


Ci 


-319.100302 


147 


Ih 


-436.702421 


185 


Ci 


-554.927247 


34 


Dsh 


-91.447319 


72 


Cs 


-205.819266 


110 


Ci 


-322.261377 


148 


Cs 


-439.807282 


186 


C2 


-558.089863 


35 


C2v 


-94.464663 


73 


Cs 


-208.794190 


111 


Ci 


-325.367188 


149 


Ci 


-442.741345 


187 


Ci 


-561.167833 


36 


Cs 


-97.324952 


74 


Cs 


-211.818219 


112 


Ci 


-328.390553 


150 


Cs 


-445.909288 


188 


Ci 


-564.330048 


37 


Cs 


-100.214276 


75 


Cs 


-215.088423 


113 


Ci 


-331.327112 


151 


C2v 


-449.086382 


189 


Ci 


-567.390087 


38 


Deh 


-103.264168 


76 


Dzh 


-218.367742 


114 


Cs 


-334.422411 


152 


Cs 


-452.190934 


190 


C2 


-570.537640 


39 


Csv 


-106.437242 


77 


Csv 


-221.514884 


115 


Ci 


-337.598496 


153 


C2v 


-455.278091 








40 




-109.401628 


78 


Dsh 


-224.595480 


116 


Cs 


-340.759024 


154 


Cs 


-458.319723 









anti-Mackay overlayer.^^ This growth maintains the poly- 
tetrahedral character of the chistcrs and leads to a se- 
ries of interpenetrating 13-atom icosahedra at iV = 19, 
23, 26, 29, 32 and 34. For example, the 34-atom struc- 
ture involves seven interpenetrating icosahedra with their 
centres in the shape of a pentagonal bipyramid. This 
growth sequence is very similar to that for Lennard- Jones 
clusters^^ although in a few instances slightly different 
sites in the anti-Mackay ovcrlayer arc occupied. 

Beyond this size is where things become different from 



any previous system. For example, for Lennard- Jones 
clusters there is a crossover at = 31 to the Mackay 
overlayer,^^ which leads to the 55-atom Mackay icosa- 
hedron. For model metal potentials, even those that 
favour Mackay icosahedra, one often sees decahedral and 
close-packed structures at sizes away from the complete 
icosahedra. Instead, for these aluminium clusters the 
polytetrahedral character of the global minima contin- 
ues. However, it is not by continuation of the anti- 
Mackay overlayer as for long-ranged Morse clusters, but 
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FIG. 4: A selection of the putative global minima. For the larger clusters the disclination network associated with a structure 
is depicted to its right. 
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through structures involving ordered arrays of dischna- 
tions. InitiaUy some of these structures are similar to that 
for a modified form of the Dzugutov potential that was 
designed to produce compact polytetrahedral clusters, 
but they then quickly diverge. 

The 40-atom global minimum is a sixfold-symmetric 
equivalent of AI34 . It involves a ring of six interpenetrat- 
ing icosahedra with two interpenetrating Z14 polyhedra 
along the symmetry axis. These latter two polyhedra give 
rise to the single disclination line that passes through the 
structure. Similar structures have previously been found 
for Dzugutov clusters, ^^'^^ and it is a small fragment of 
the Frank-Kasper Z phase. 

Similar clusters are seen for = 44 and 48, but with 
a Z15 and Z16 polyhedron, respectively, lying slightly 
off-centre; these act as three-fold and four-fold nodes in 
the disclination networks. Only for AI48 does further 
growth occur on this cluster leading to the particularly 
stable 61-atom global minimum with tetrahedral sym- 
metry. This cluster was first identified for a long-ranged 
Morse cluster, -'^'^ and has since been found for clusters 
with a generalized set of metallic interactions,'^* and in 
experiments on strontium clusters.* 

Two exceptions to this dominance of polytetrahedral 
structures occur at = 54 and 55, where a complete 
Mackay icosahedron and one with a vertex missing atom 
are most stable. However, away from this magic number 
structures based on Mackay icosahedra quickly become 
significantly higher in energy than the polytetrahedral 
structures. 

The 64-atom global minimum is the first to involve 
not just a single network of connected disclinations. It 
is based on Algi but it also has a single linear disclina- 
tion passing through the clusters. Such disclination lines 
are a common motif for the larger clusters. At A^ = 69 
the global minimum has an unusual structure with a 19- 
coordinate central atom that acts as a node for six discli- 
nations going out from it. 

Starting at N—72 a new class of polytetrahedral clus- 
ters becomes most stable. They are all based in some 
way on a set of Frank-Kasper phases, which are layer 
structures'^ and can be considered as a tiling of squares 
and triangles that are suitably decorated with atoms. 
The four simplest crystalline tilings are represented in 
Figure 5(a). The four Frank-Kasper phases correspond 
to the four possible ways in which a point can be sur- 
rounded by squares and triangles. The Z phase has six 
triangles around each vertex, the A15 phase four squares 
and the H and a phases correspond to the two ways that 
two squares and three triangles can share a common ver- 
tex. 

Overlaid on the tilings in Figure 5(a) are the associated 
disclination networks. All have planar arrays of disclina- 
tions that are mediated by Z14 and Z15 polyhedra and 
that are threaded by linear disclinations that are per- 
pendicular to the plane of the tiling. A number of the 
global minima can be understood simply in terms of finite 
square-triangle tilings, as illustrated in Figure 5(b). Al7g 



(a) Z phase Al 5 phase 




FIG. 5: (a) The four square-triangle Frank-Kasper phases. 
The diagrams show a portion of the square-triangle tiling 
(thin black lines) for these phases, (b) Some of the square- 
triangle motifs that occur in the clusters in Fig. 4 labelled 
by their sizes. The disclination networks (thick red lines and 
dots) associated with these tilings are also depicted. The 
red dots represent disclinations coming out of the page. The 
disclination lines that cross in the square tiles are at different 
heights in the crystal. 



has a disclination network that is part of the Z phase; 
AI72 and Alga correspond to the H phase, and N = 138, 
142, 160 and 172 to a cr-like tihng. The latter slightly de- 
viates from the a phase, because, although the local ar- 
rangement of squares and triangles around all the vertices 
in the tiling is the same as for the a phase, the central 
triangle only shares edges with squares. In fact this motif 
is common in models of dedecagonal quasicrystals.^*"^° 

The other larger global minima can be basically 
thought of in terms of these tilings, but with some devi- 
ations in the disclination networks at the cluster surface. 
So A = 107, 120, 132, 134 and 148 have H-like dischna- 
tion networks, where the Algs tiling has been extended 
downward for A^ = 134 and across for A^ = 120, 132 and 
148. Similarly, A^ = 94, 116, 151 and 186 also have cr-like 
tilings. All the a-like structures are heading towards or 
based upon Aligo. 

The deviations from the square-triangle structures in- 
volve breaking the planarity of the disclination networks 
to connect disclinations in different layers of the struc- 
ture. The most common is to link two parallel disclina- 
tions that would otherwise have simply exited the clus- 
ter. This connection gives a disclination pattern similar 
to that for a three triangle tiling, but in a plane perpen- 
dicular to the usual plane of the disclination network. 
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For example, the only difference between the structures 
at TV = 94 and 116, and 138 and 151 is the addition of 
one of these linkages. Similarly, by the addition of two of 
these linkages one can generate Alige from Align- This 
motif is also present in the structures at iV = 120, 132, 
134 and 148. The other linkage is similar and involves 
linking what would otherwise have been three-fold nodes 
in the disclination network, thus generating Z16 nodes. 
This type of connection occurs for A^=94 and 116. 

It is noticeable that the variations in the energy with 
respect to the average for clusters with more than 100 
atoms are relatively small (Fig. 3). This is because there 
are lots of ways of creating polytetrahedral packings. 
This behaviours is unlike what occurs, say, for struc- 
tures based on the Mackay icosahedra, where, between 
the magic numbers, structures with incomplete outer lay- 
ers give rise to a substantial energy variation. Even the 
most stable large polytetrahedral packings, the cr-like 
structures centred around N — 166, only give rise to 
a broad minimum in Fig. 3. The structure at the mini- 
mum, Aligg, is in fact intermediate between the 160- and 
172- atom structures depicted in Fig. 4 with three of the 
six disclinations extended. 



IV. ALei: A CASE STUDY 

To try to unravel why the aluminium glue potential 
favours polytetrahedral structures, in this section I take 
Algi as a case study and compare the energetics of some 
competing structural forms (Table II). The structures 
I consider are the polytetrahedral global minimum (Fig. 
4) , an icosahedral structure that can be formed by adding 
a six-atom overlayer to the 55-atom Mackay icosahedron 
and an fee cluster. The latter two are illustrated in the 
insets to Fig. 6. 

The decomposition of the potential energy into glue 
and pair contributions shows (Table II) that the global 
minimum is stabilized by a significantly lower .Epair- This 
is partially offset by a higher -Egiuc than the icosahe- 
dral structure. Figures 6(a) and (b) enable us to un- 
derstand why the pair energy is lowest for the polytetra- 
hedral clusters, where (a) depicts npairs(?'), the number 
of pairs of atoms that are separated by less than r, and 
(b) impair (j'ij < f) , the contribution to E'pair from pairs of 
atoms that are separated by less than r. 

The contribution to impair from the nearest neighbour 
shell (defined as those pairs with a separation less than 
3.5A) can be decomposed into two terms 

i?p\"i, = nn„0((r„„)) + A£;pt', (4) 

where rinn is the number of nearest neighbours, (rnn) is 
the average separation of pairs in the nearest-neighbour 
shell. The first term is the impair value that would re- 
sult if all the rinn nearest neighbours had the same pair 
separation, and the second term, AE'p^'j*, is the energy 
cost associated with the nearest-neighbour pair distances 
having a distribution of values. 
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FIG. 6: A comparison of the properties of the global min- 
imum (gmin), the lowest-energy icosahedral structure (icos) 
and the lowest-energy fee structure for N — 61. The latter 
two structures are illustrated in the corners of (b) and (c) , re- 
spectively, (a) npairs(r). (b) Ep^^ir^Vij < r). (c) {p{rij < r)). 
In (a) vertical lines have been added corresponding to the 
values of r where the pair potential passes through zero, is a 
minimum or a maximum, and goes to zero at the cutoff, as 
labelled. 
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TABLE II: The properties of a series of 61-atom structures, 
namely the global minimum (gmin), the lowest-energy icosar 
hedral (icos) and lowest-energy fee structure, ffp is the stan- 
dard deviation in the values of pi, and ct,. is the standard de- 
viation in the nearest-neighbour distances, n-buik is the num- 
ber of atoms in the interior of the cluster, where an atom is 
deemed to be on the surface if its coordination number is ten 
or fewer. -Bbuik and -Esurf are the average energies of the atoms 
in the interior and on the surface of the cluster, respectively. 
All other quantities are defined in the text. 





gmin 


icos 


fee 


point group 


Td 






E / cV 




-171.697 


-170.975 


-Bpair / eV 


-22.795 


-20.901 


-20.991 


i^giue / eV 


-149.992 


-150.796 


-149.984 


Ell,. 1 eV 


-5.899 


-5.284 


-5.680 


^nn 


284 


260 


249 


(rnn) 


2.8536 


2.7973 


2.7974 


(7r 


0.1565 


0.0742 


0.0442 


AS^if, / eV 


4.817 


0.997 


0.342 


(P> 


0.7078 


0.7098 


0.6906 




0.2145 


0.2035 


0.1904 


AE^iil 1 eV 


6.052 


5.338 


5.159 


n-bulk 


17 


14 


13 


^'buik / eV atom~^ 


-3.104 


-3.119 


-3.186 


^-surf / eV atom~^ 


-2.728 


-2.724 


-2.699 



As is well-known icosahedral clusters have larger val- 
ues of film than fee elustcrs, because Mackay icosahedra 
are quasispherical and only have high-coordinate {111} 
faces, whereas the most spherical fee clusters must have 
a significant proportion of lower-coordinated {100} faces. 
This more favourable surface energy is the reason why 
icosahedral clusters are often lowest in energy for small 
clusters interacting with pair potentials. The polytetra- 
hedral global minimum has an even higher value of rinn 
than the icosahedral structure because of the higher coor- 
dination number for the interior atoms that lie on discli- 
nations, and the densely packed surface. 

As anticipated in Section II A, (rnn) is significantly 
shorter than the distance corresponding to the first mini- 
mum in the pair potential — some nearest-neighbour pairs 
in the polytetrahedral and icosahedral structures even 
have positive pair energies, (rnn) also has a significant 
structural dependence. For the icosahedral and fee struc- 
tures, (rnn) is contracted with respect to rnJf^', as is usual 
for metal clusters, in order to try to increase the p val- 
ues of the surface atoms. By contrast, (rnn) for the global 
minimum is actually slightly larger than rnJf'. This is not 
that there is no contraction at the surface, but rather that 
this contraction is with respect to the appropriate bulk 
Frank-Kasper phase. 

Given the large value of rinn for the global minimum 
and that (rnn) lies closest to the potential minimum, it is 
not so surprising that the global minimum has the lowest 



value of -Epfii-. However, it is only marginally lower than 
that for the fee cluster, which in turn is lower than the 
icosahedral cluster (Figure 6(b)) even though it has fewer 
nearest neighbours. These effects can be understood by 
taking account of the distribution of values. As (/)(r) 
has positive curvature in this region, AiJp^'?* must be pos- 
itive. 

As is evident from Figure 6(a) and the values in 
Table II, the global minimum has the widest distribu- 
tion of nearest-neighbour distances, and hence by far the 
largest Ai^p^'-*. By contrast the fee cluster has the nar- 
rowest distribution and smallest AiJpif;*. This explains 
why £^pair for the fee cluster is actually lower than that 
for the icosahedral cluster and only just above that for 
the global minimum. The large value of for the global 
minimum reflects the distortions of the tetrahedra away 
from regularity that are necessary in any polytetrahedral 
packing. 

E'p^i^, however, only represents roughly a quarter of 
the total pair energy; the dominant contribution energy 
comes from next neighbours. There are a number of 
possible configurations for such neighbours. The near- 
est next neighbours correspond to those at opposite ver- 
tices of an octahedron, which if regular corresponds to 
a separation of a/2 rnn = 1.4142 rnn- These configura- 
tions occur most frequently in the fee clusters, also occur 
in the twenty fee tetrahedra that make up the Mackay 
icosahedron, but do not occur at all in a polytetrahedral 
structure. 

The next set of pair distances correspond to atoms 
at opposite ends of a trigonal bipyramid, wh ich i f it is 
made of two regular tetrahedra occurs at 2-\/2/3rnn = 
1.6330 rnn- These configurations are absent from fee 
clusters, occur around the fivefold axes of the Mackay 
icosahedra and also across the twin planes between the 
strained fee tetrahedra that make up the icosahedron and 
are most prevalent in polytetrahedral clusters. 

Then, there are pairs at the extreme ends of a planar 
rhombus, which if it consists of two equilateral triangles 
would be at a separation of \/3rn„ = 1.7321 rnn- How- 
ever, these pairs, which are most common in the fee struc- 
tures due to the close-packed planes, are less important 
to the energetics because, of the three sets of configura- 
tions, they are the furthest from the second minimum in 
the pair potential. 

Given that rniin2 = 1.509 rnJf', one might expect from 
the above ideal configurations that the octahedral pairs 
would lie slightly closer to the second potential minimum 
than those pairs at opposite ends of a trigonal bipyra- 
mid. However, the contraction of (rnn) with respect to 
the crystal takes the octahedral distances further away 
from this minimum. Furthermore, the wide distribu- 
tion of nearest neighbour distances in the polytetrahedral 
structures also leads to a wide variation in the dimensions 
of the trigonal bipyramids. As can be seen from Figure 
6(a) a significant fraction of these polytetrahedral pairs 
for the global minimum lie closer to the pair potential 
minimum than the octahedral distances. These effects 
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are responsible for the substantially lower pair energy of 
the global minimmn (Figure 6(b)). The trigonal bipyra- 
mids present in the ieosahedral structure also cause the 
difference in the total pair energy between it and the fee 
structure to narrow, compared to that for -Ep^ir- 

Similar to the analysis for E^^ir the glue energy can be 
written as 

E^,,, = NUm + AEX. (5) 

where the first term is the glue energy that would result 
if all atoms had the same value of p {{p) is the average 
value of p), and the second term is the correction due to 
the {pi} having a distribution of values. 

We can examine the role of different pair distances on 
the glue energy by following {p{rij < r)), the contribu- 
tion to (p) from pairs of atoms that are separated by 
less than r. It is clear from Figure 6(c) that, as antic- 
ipated, the major contribution to {p) comes from near- 
est neighbours, and the values of {p{rij < r)) after the 
nearest-neighbour shell reflect the relative values of Unn- 
(rnn) is also important. The average contribution to (p) 
from each nearest neighbour in the global minimum is less 
because of the longer (rnn)- Indeed, if the three struc- 
tures had the same (rnn) the polytetrahedral structure 
would have comfortably the lowest glue energy. The dif- 
ferences in CTr have little effect, because the function p(r) 
is approximately linear for the relevant range of r (Figure 
2(a)). 

As p{r) is a monotonically decreasing function, the 
next largest contribution to (p) comes from octahedral 
next neighbours. All other pairs have little effect on the 
relative values of (p) for the different structures because 
by these distances p{r) is small in magnitude. 

The effect of next-nearest neighbours is to cause (p) 
for the ieosahedral structure to just overtake that for the 
global minimum, and similarly to diminish the difference 
between (p) for the global minimum and the fee structure, 
compared to the relative values of (p(rjj < r)) after the 
nearest-neighbour shell. 

The ordering of the (p) values agrees with that for 
£^giue. However, although the global minimum has a sim- 
ilar value of {p) to that for the ieosahedral structure, its 
glue energy is only just greater than that for the fee struc- 
ture. This effect can be explained by taking into account 
the distribution of pi values. As U{p) has positive cur- 
vature, Ai?gj^^* must be positive, and is in fact largest 
for the global minimum because it has the greatest vari- 
ance in Pi values (Table II) . To a first approximation the 
distribution of pi values is bimodal with similar values 
for surface atoms, and similar values for interior atoms. 
As the polytetrahedral cluster has the most atoms in the 
interior of the cluster due to the high coordination num- 
bers, it has the largest ap. 

As we have pointed out above the value of (rnn) is 
important for the values of both impair and E'giuei so it 
is important to understand what determines (rnn). At 
(rnn) the decrease in the glue energy by uniformly shrink- 
ing the cluster is exactly balanced by the increase in 



the pair energy. This energy balance occurs at larger 
(rnn) for the global minimum for a number of reasons. 
Firstly, the decrease in the glue energy as the cluster is 
shrunk (this increases the p values of the surface atoms) 
is smallest for the global minimum because it has the 
fewest number of atoms in the surface of the cluster and 
has some atoms with high coordination number in the 
cluster interior (see below). Secondly, the non-linear 
increase in the pair energy as the cluster is shrunk is 
largest for the global minimum. This mainly stems from 
the nearest-neighbour contribution and reflects both the 
greater number of nearest neighbours for the global min- 
imum and its greater ar- 

The average nearest-neighbour distance also shows a 
systematic dependence on the coordination number. For 
energetic (see Section II A) and geometric reasons the 
atoms with higher coordination numbers have larger av- 
erage nearest-neighbour distances. For example, the 
nearest neighbours of the sixteen-coordinate atom in the 
centre of the global minimum are at an average separa- 
tion of 2.933 A, whereas the corresponding distance for 
the 12-coordinate atoms is 2.725 A. 

Although the analysis in this section has been done 
here for a single case study, similar analyses for other 
cluster sizes have revealed the same principles at work. 



V. LARGER CLUSTERS 

Although polytetrahedral clusters are dominant in the 
size range considered in the last section, one would expect 
that clusters at sufficiently large size would adopt the 
bulk structure, which for aluminium is fee. This assumes 
that the current glue potential can correctly reproduce 
the bulk structure. However, to the best of my knowl- 
edge, the relative energetics for different crystal struc- 
tures have not been compared for this potential; for a 
somewhat similar aluminium empirical potential the fee 
crystal is only marginally lower in energy than an A15 
Frank-Kasper phase. ''^ 

That for Algi, and all the other clusters we examined 
in detail, £^buik lower than that for the polytetrahe- 
dral and ieosahedral clusters (Table II) suggests that the 
system will converge to the correct bulk structure. To 
further understand the size evolution of the cluster struc- 
ture I have reoptimized a series of Mackay icosahedra, 
fee truncated octahedra and Marks decahedra, which are 
usually the optimal shapes for these morphologies.^^ The 
energies of these sequences are compared in Figure 7(a) 
along with the global minima obtained in the Section III. 
I have not attempted to construct a sequence of larger 
polytetrahedral clusters, because one cannot easily de- 
duce from the global minima how the polytetrahedral 
growth sequence should continue. Instead, we also plot 
i?ave, the fit to the energies of the global minima. 

From Figure 7(a) it is clear that decahedral structures 
are never competitive for this potential. This is slightly 
unusual as it is often found that for materials which are 
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close-packed in bulk, there is a progression from icosahe- 
dral to decahedral to close-packed structures as the size 
increases. Instead, the icosahedral and fee lines cross at 
N PS 520 suggesting that beyound this size fee clusters 
will be lower in energy for the majority of sizes. The i^ave 
curve also seems to suggest that these fee clusters will 
be global minima, but one has to be somewhat careful. 
Firstly, one should not necessarily trust this function out- 
side the size range to which it has been fitted. Secondly, 
it represents the average energy of the polytetrahedral 
global minima, whereas the other curves in Figure 7(a) 
are lower bounds to the energies of the clusters of that 
type. For example, £^ave only dips below the interpolated 
energies of the complete Mackay icosahedra for a small 
range of A'', even though polytetrahedral structures are 
the global minima for the vast majority of the size range 
considered in Section III. 

To get a better idea of the causes of this size evolution 
I have plotted in Figures 7(b) and (c) the glue and pair 
contribution to the energies of these sequences of struc- 
tures. It is immediately clear, in agreement with the 
analysis of Algi, that the polytetrahedral global minima 
arc stabilized by their particularly low pair energies. Fur- 
thermore, the global minima that are based on Mackay 
icosahedra at A/' = 54-55 and 147 clearly stand out from 
this trend. 

It is also clear that the crossover from icosahedral to fee 
structures is driven by the increasing advantage provided 
by its lower pair energy. By contrast, the glue energy of 
the icosahedron only becomes greater than that for the 
truncated octahedra at A'' « 6500. E^^^ more strongly 
reflects the greater number of nearest neighbours asso- 
ciated with the Mackay icosahedra and is less sensitive 
to the strain inherent to the icosahedra. By contrast, as 
noted in section III and as is usual for pair potentials,^'' 
the nearest neighbour contribution to the pair energy is 
reduced due to the strain in the pair distances for the 
icosahedra. 



VI. EXPERIMENTAL AND THEORETICAL 
COMPARISONS 



-3^ . . . 

10 100 ^ 1000 10000 

FIG. 7: A comparison of the relative values of (a) E, (b) 
impair and (c) iSgiue for the global minima for < 190 
(gmin) , aud sequences of fee truncated octahedra with hexag- 
onal {111} facets (fee), Mackay icosahedra (icos) and Marks 
decahedra (deca). In each graph the energies are mea- 
sured relative to fits to those for the truncated octahedra, 
where Btoct = -3.360iV -|- LSgOAT^/^ -f X.Vjm^l^ - 0.267, 
-Bfoct = -2.657iV + 0.456iV^''^ 1.096iV^''^ + 0.621 and 
Elf^ = -0.703A^ -I- 1.435Ar2/3 + omSN^^^ - 0.873. In (a) 
Save, the fit to energies of the global minima used in Fig. 3, 
has also been plotted. 



Although I have focussed on the general principles that 
can be drawn from the current results, given that the 
potential has been successfully used in a wide variety 
of applications, with in some cases excellent agreement 
with experimental data,^^ one would hope that this sys- 
tem would also be a reasonable model for real aluminium 
clusters. However, one should remember that to cor- 
rectly predict the energetics of cluster structure requires 
that the potential can successfully describe a wide variety 
of surface and bulk properties for a variety of structural 
types, and so is an extremely stringent test of a potential. 

Although aluminium clusters have been much stud- 
ied experimentally, there is little reliable information 
about the geometric structure in the size range of the 
global minima considered in Section III. The difficulty 
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is that features in the experimental data can originate 
from the electronic shell structure rather than the ge- 
ometric structure. Indeed, following the successful in- 
terpretation of many of the properties of alkali metal 
clusters in terms of the electronic shell structure, alu- 
minium, a nearly free-electron metal, was a natural sys- 
tem to see how widely applicable these ideas might be. 
Indications of this electronic shell structure were found 
from the ionization potentials,®''* the magic numbers in 
mass spectra,®^'®** cluster ion mobilities®^ and photoelec- 
tron spectroscopy®®"®^. 

However, a significant number of the features in these 
properties could not be explained thus, and so these 
were usually attributed to the effects of the geometric 
structure. Noteworthy of these is the magic number at 
iV = 55, whose most obvious interpretation is in terms of 
a Mackay icosahedron. No convincing explanations, how- 
ever, have been given for the other 'geometric' features. 
For example, except for N = 55, the magic numbers do 
not correspond to those of any of the usual structural 
forms for clusters, but neither do they correspond to the 
new polytetrahedral magic numbers found here. 

Theoretically, aluminium clusters in this size range 
have been recently studied using two empirical potentials, 
namely the Murrell-Mottram®^'™ and Gupta'^^ forms. 
For both these system icosalicdral structures predomi- 
nate for the global minima, although with fee and decahe- 
dral structures competing at sizes intermediate between 
the complete Mackay icosahcdra. I have also rcoptimized 
a series of icosahcdral, decahedral, fee and polytetrahe- 
dral clusters for the potential of Mishin et ai, whose 
(j){r), p{r) and U{p) functions resemble those for the cur- 
rent glue potential and were fit to a wide range of bulk 
data.^^ In this case polytetrahedral structures again pre- 
dominate. 

There have also been many electronic structure cal- 
culations of small aluminium clusters using a variety 
of methods. ®®'®*'^^"^* For small clusters there has been 
some success in identifying the structure for A'' < 
20 by comparing calculated photoelectron spectra with 
experiment.®®'®® However, for somewhat larger clusters 
there is little consensus among the results, partially be- 
cause the computational expense prevents one rigorously 
searching for the global minimum — usually only a selec- 
tion of clusters are compared — but also because different 
methodologies simply give different results. For example, 
there is disagreement over whether for AI55 the Mackay 
icosahedron^® or the fee cuboctahedron®®'^"^'^^ is lower in 
energy, and still others have found disordered structures 
to be lower in energy than both.'^® A similar range of 
results have been found for Ali47.'''^'^^ 

For large aluminium clusters there are clear experimen- 
tal indicators of the geometric structure. For N > 250 
(up to at least N « 10 000) there arc set of regular 
spaced magic numbers in the mass spectra.®**'^^"®^ Ini- 
tially, this result was interpreted in terms of electronic- 
shell effects, ®^'®^ but Martin et al. soon showed that the 
magic numbers were in fact due to octahedra, with each 



additional magic number corresponding to the addition of 
a complete overlayer to a single face of the octahedra. ^^'®^ 
Growth simulations further added to the plausibility of 
this interpretation.®'*"®''' 

For large clusters one should always ask whether the 
observed structures occur because they are most stable 
or because of the growth kinetics. This is especially the 
case here, because the most stable fee clusters are usu- 
ally more spherical truncated octahedra.®® rather than 
octahedra. Indeed as the life times of the clusters are in- 
creased, the magic numbers disappear,®*^ indicating their 
kinetic origin. However, it is probably reasonable to con- 
clude that fee structures are most stable in this size range, 
even if the appearance of octahedra is a kinetic effect that 
results from the miclcation of a new face being slower 
than the subsequent growth. 

My estimate that fee structures become lowest in en- 
ergy beyond A'' k, 500 is therefore reasonably close to ex- 
periment. This contrasts with the Gupta potential where 
fee structures only become more stable beyond N w 1250 
and the Murrell-Mottram potential where the crossover 
occurs at even larger sizes.®'' However, one should remem- 
ber that these energetic crossover sizes are only strictly 
valid at T = OK. If one compares the relative free en- 
ergies, the crossover sizes can be sensitively dependent 
on the temperature, mainly due to the differences in the 
vibrational entropy of the competing structures.®®'®^ 



VII. CONCLUSION 

I have found that the global minimum of aluminium 
clusters described by a glue potential to be predomi- 
nantly polytetrahedral for N < 160. The structure of the 
larger clusters can be described in terms of their discli- 
nation networks, which resemble those of the square- 
triangle Frank-Kaspcr phases. This result is somewhat 
surprising given that ordered polytetrahedral structures 
are only found for alloys in bulk, but it reflects the greater 
tendency for polytetrahedral order in clusters. The bulk 
behaviour should of course emerge for sufficiently large 
clusters, and our results suggest that for N > 500 clusters 
with the bulk fee structure predominate. 

An important aspect of this work has been to relate 
the observed structures back to the interatomic interac- 
tions. Polytetrahedral structures arc favoured for a com- 
bination of reasons. Firstly, the clusters present surfaces 
with a low surface energy and high average cooordination 
number, so will be potentially favoured at small sizes. 
Secondly, the glue function, which gives rise to most of 
the cohesive energy of the clusters, is relatively insen- 
sitive to the strain in the pair distances and the high 
coordination numbers that are inevitable in polytetrahe- 
dral structures. Thirdly, as with the modified Dzugutov 
potential,^® the two minima in the pair potential most 
closely match the distances in polytetrahedral structures. 

As the first two of these reasons arc true for most metal 
potentials, it is not so surprising that there is experimen- 
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tal evidence of polytetrahedral structures for monatomic 
cobalt, strontium,^ and perhaps some rare-earth^^ clus- 
ters. Nor would it be surprising if the number of theo- 
retical and experimental examples were to increase. My 
results further^^ illustrate the types of polytetrahedral 
structures that one might expect for such clusters, along 
with a series of magic numbers that will be useful in com- 
parisons with experiment. 

Whether the preference for polytetrahedral structures 
in the current model carries over to real aluminium clus- 
ters is not yet clear. Because of the lack of easily inter- 
prctablc experimental data on the geometric structure of 
aluminium clusters with less than 190 atoms, it is hard to 
draw any firm conclusions about the realism of the struc- 
tures predicted by the glue potential in this size range. 
Moreover, ab intio calculations add little further insight 
because of the lack of consensus. However, comparisons 



to the experimental data for larger clusters suggest that 
the glue potential is more realistic than the other em- 
pirical potentials for which calculations have been per- 
formed. 

The properties of these clusters can also provide use- 
ful insights into the behaviour observed for this poten- 
tial in other contexts. For example, studies on alu- 
minium nanowires have found polytetrahedral structures 
for those wires with sufficiently small cross-section.*^ 
Some of these involve only icosahedral local coordina- 
tion, but others clearly involve disclinations. For exam- 
ple, the wire that Giilseren et al. denote as A7 has a 
single disclination running along the centre of the wire. 
Similarly, from the current results it should be unsurpris- 
ing that there is strong polytetrahedral ordering in the 
aluminium melt described by this potential.*^ 
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